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Abstract 

Using strong lensing data Milgrom's MOdified Newtonian Dynamics (MOND) or its covariant 
TeVeS (Tensor- Vector-Scalar Theory) is being examined here as an alternative to the conventional 
ACDM paradigm. We examine 10 double-image gravitational lensing systems, in which the lens 
masses have been estimated by stellar population synthesis models. While mild deviations exist, 
we do not find out that strong cases for outliers to the TeVeS theory. 
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I. INTRODUCTION 



The concordant ACDM paradigm [see e.g., Ill has been accepted as a successful framework 
to understand the evolution of the universe. However, in this standard model there are two 
unexplained dark sectors, dubbed dark matter and dark energy, still need more fundamental 
understandings. One of various endeavors on this explanation is the approach of modified 



gravity, including Bekenstein's TeVeS [2]. Unlike most theories of modified gravity that 
are proposed for explaining dark energy, TeVeS was originally proposed to explain dark 
matter instead. Indeed, TeVeS was built up by Bekenstein for a viable relativistic version of 
Milgrom's MOdified Newtonian Dynamics (MOND) [3]. Other many recent theories which 
recovers MOND in it's limit include the bimetric theory of {4] and the dark fluid theory 

of Q. 

In MOND, the demand of the exotic dark matter is replaced by the modification of 
Newton's second law: 

/2(|a|/ao)a= -V$ N , (1) 

with do ~ 1.2 x 10 10 m s~ 2 . p,(x) is called interpolation function in literature sometimes. 
jx(x) ~ x for 1 C 1 (the deep MOND regime), and fx{x) — > 1 when x — > 1 (the Newtonian 
regime). Here x = |a|/oo, the ratio of the acceleration to do, is a measure of modified gravity. 
MOND -in its relativistic version TeVeS |2l-is not only able to explain CMB successfully 

n q 

with 2 eV neutrino mass [6|, |7|] but also even more successful than the ACDM paradigm on 
the dynamics of spiral galaxies jsl, 9]. 

Unlike cold dark matter, however, massive neutrino cannot aggregate at galactic scales 
and below, so any evidence for the demand of dark matter at these scales is devastating 
to MOND. Beside dynamical analysis on rotation curve (of spiral galaxies), gravitational 
lensing offers us another way to investigate MOND and TeVeS at these scales. Indeed, 
the requirement of dark matter in the bullet cluster makes gravitational lensing even more 
important flo| |. 



The earliest works about gravitational lensing in MOND were studied by [ll|] and [12 ]. 
Although their works were performed even before the appearance of a viable general covari- 
ant MOND, and their calculations on the angle of deflection were artificially forced into the 
deep MOND regime when the acceleration is less than Oq, i.e., fl(x) = x for x < 1 (not 



x <§C 1), [ll| and [12| could show a good guess on how to calculate lensing in MOND. After 



the advent of TeVeS, a first calculation on light bending in TeVeS was done by Bekenstein 



himself, and was followed up by 13]. Based on a point mass model, [l3| apply the for- 
malism of gravitational lesing in TeVeS to the theoretical discussion on angle of deflection, 
magnification, as well as time delay. 

On the other hand, taking from a more phenomenal approach, [l^] showed that TeVeS is 
consistent with a sample of double-image lenses in CASTLES catalogue by modeling lenses 



as Hernquist model. 15[ studied the effects of asymmetric systems on gravitational lensing. 
li| then applied non-spherical model to investigate strong lensing in TeVeS, and found 10 
out of 15 systems are consistent with TeVeS. All of other 5 systems are found to reside in or 
close to clusters of galaxy, where external field could have significant influence. Furthermore, 
the effect of large filaments on gravitational lensing was studied by [17]. They argued that 
filamentary structures might have complex but significant contribution to the system such 
as bullet cluster, so the need of dark matter in bullet cluster might be spared again in 
MOND. However, all of these work above are non-relativistic approximation of TeVeS. In a 
conference paper, [3] also showed the lens data from CASTLES and SLACS catalog are 
consistent with TeVeS, but they did not show the details. 



The first effort to calculate gravitational lesing from first principle was given by 32]. 



Along with their other works 



19 



20] , they showed that TeVeS might still need dark 



matter to explain the lensing systems from CASTLES, and is lack of consistent results 
between dynamical and gravitational lensing analysis fiol . Q] . 

Due to the importance of studying gravitational lesning in TeVeS, and the inconsistency 
of the results on this field, in this paper, we are goin g to apply the relativistic approach 



developed by 02Bekenstein [2j and 04Chiu et al. [13 1, and the Hernquist model (used for 
TeVeS firstly by 34Zhao et al. [ijjjto the 10 double-image systems from a total of 18 
objects studied in 12Ferreras et al. [21]. In this paper we will also clarify the underlying 
assumptions of lensing calculation in TeVeS. The structure of the paper is organized as 
below. In Sec. [Ill we will briefly outline the formalism of gravitational lensing in TeVeS. 
We will also discuss how gravitaional lensing in TeVeS will depend on different choice of 
/}(|a|/ao), and its application to double-image systems. We then will give our result in 
Sec. IHIl and discussion in Sec. [TV] 



II. FORMALISM OF GRAVITATIONAL LENSING IN TEVES 



A. Lensing in TeVeS 



While considering strong lensing in TeVeS, we take the assumption that in the weak field 
limit, the physical metric of a static spherically symmetric system can be expressed in the 
isotropic form (c=l): 



9o 



■ dx a dx 13 



[1 + 2<$>)dt 2 + (1 - 2§)[dg 2 + Q 2 {d6 2 + sin 2 



0] 



(2) 



where $ = S$jy + 4> and H = e _2 ^ c (l + K/2)~ l with <ft c as the asymptotic boundary value 
of 4> and K < 1CT 3 from the constrain of PPN parameters ^]. Since in TeVeS all kinds of 
matter are coupled to physical metric rather than Einstein's metric, the connection in the 
geodesic equation, 

(3) 



d 2 x^ dx v dx x 

+ r^- — — = o 



dp 2 ' uX dp dp 

has to be constructed from physical metric. Here p is some proper affine parameter. In 
the case of light bending by a static spherically symmetric lens, we can apply Eq.(j2]) to the 
geodesic equation, Eq.()3]). Consider a photon propagates on a null geodesic, and moves in 
the equatorial plane, i.e., 9 = ir/2, the geodesic equation or equations of motion can be 
written as 



1 + 2$)i = E 



(1-2$)$V 

d - 2$)q 2 + (1 + 2§)q- 2 L 2 - (1 - 2<&)E 2 







(4) 
(5) 
(6) 



where over dot denotes the derivative with respect to p. Recall the fact that at the closest 
approach g = 0, q = Qo, we can express the ratio of the angular momentum L to energy E 
as 

2 = L 2 gg(l - 2$ ) 
" E 2 (l + 2$ ) ' 1 } 

where $ = ^(f?o)- Here b is called the impact parameter. Combining Eqs. ©-© gives the 

orbit of the photon, 



1-4$) + (1 -4$ c 



Qo 



dg 
g 2 \d(p 



+ 1 



0. 



(8) 
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of which the solution in quadrature is 




<p= I <[£) a [l-4(*-$o)]-l| j- (9) 

If we take the Taylor expansion of Eq. to the first order of $, we get 

^ 2g dg + ^ r 4g($-$ ) ^ 

The first term is the orbit without gravity, i.e., a straight line (<p — > it as g — > oo). The 
second term is the angle of the deviation due to the gravity. Hence, to first order of <3>, the 
deflection angle is 

M^($-$ ) f s 4|V$| , 4^o$ 



-1/2 



where we have made use of integration by parts. Since for very large g, $ behaves as In g, it 
is legitimate to ignore the second term of Eq. lfTTj) for strong lensing systems [3]. We then 
have the kernel equation of strong gravitational lensing in TeVeS (up to first order of V$), 

r°° |v$| 

w^ ie - (12) 

Recall that go is the distance of the closest approach, and V$ is the MONDian gravity. 



B. Gravity in TeVeS 

In last subsection we have shown that, under the assumption of the validity of Eq.Q, 
the difference between the angle of deflection in GR and that in TeVeS only arises from V$. 
In this subsection, we are going to discuss the relation between Newtonian gravity, V<3?at, 
and MONDian gravity, V$, and its application to strong lensing. 

It has been shown in quasi-static limit, the MONDian potential in TeVeS can be expressed 
as the combination of Newtonian potential and a scalar field 

$ = S$ w + 0, (13) 

where S is a parameter in TeVeS and is approximately 1, and <fi represents the strength of a 
scalar field. Moreover, the scalar field itself is linked to the Newtonian potential via a free 
function fi [2j, 

V</>=(k/4irfj,)V$ N . (14) 



where /x is a function of |V0|. This free function \x should be chosen carefully in order to 
reproduce Newtonian or MOND behavior at quasi-static limits. In fact, Eq. flT3|) . Eq. ffT4j) . 
and Eq. flTJ give the relation between \i and p, of the modified Poisson equation or Milgraom's 
law Eq. flTJ) 

(15) 



1 _ „ k 
jl Anfi 



TeVeS has only two parameters (k and 5) and one free function (//). Therefore, the MON- 
Dian behavior controlled by jl in Milgrom's law Eq. ([1]) could be understood via the param- 
eters and free function in TeVeS. On the other hand, we are able to express the modified 
gravity V$ (even in TeVeS) as a function of Newtonian gravity V$at via jl. In the following, 
we are going to discuss three commonly used interpolation functions, jl. 

In Bekenstein's TeVeS paper he proposed the following interpolation function 

-l + yr+4i 



jl{x) 



(16) 



where x 
galaxies 



1 + VI + 4x ' 

a|/oo. The Bekenstein's form above fails to fit the rotation curves of spiral 
22J. In order to apply the MONDian lens equations to observational data, we 
also study the simple form and the standard form which fit the rotation curves better. The 
simple form is [3] 

x 

f l ( x ) = T—:' ( 17 ) 

and the standard form is 

jl(x)= j=^== . (18) 



1 + x 



X 



Even though the standard form will lead to bi- values problem, so is supposed to be unphysical 
under the framework of TeVeS 



23 



24| , here we still compare the standard form with the 



other two by treating it as empirical function and for the sake of curiosity. 

In fact, all these three forms can be included in the following two-parameter form 

jl(x) 



(19) 



(1 + r]x a ) + a/(1 - r)x a ) 2 + 4x c 

Here, (a, rj) = (1,0), (1,1), (2,1) and (oo, 1) correspond to Bekenstein form, simple form, 
standard form and the naive sharp-break form jl = min(l,x), respectively. Recall that we 
have set c = 1. We note that 31 Zhao & Famaey [24] have proposed a similar expression in 
which they combined Bekenstein form and simple form. We may call Eq. ( !T9|) the invert- 
ible canonical interpolation function which goes from the naive sharp-break to the smooth 



6 



Bekenstein form. The corresponding /x can be found by Eq. (|T5|) . The nicest thing of our 
invertible canonical interpolation function is that it has a very simple counterpart in the re- 
cent Quasi-MOND theory or its relativistic version called Bi- metric MOND {4, 
Eq. (jTJ) with jx given by Eq. ffTO]) gives 



251 ] . Inverting 



— a 



iv 



(20) 



where 



u(x N ) 



l/a 



(21) 



and xn = |V$at|/oo- This ana! 



Note that some earlier 
law as 
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ytical result can greatly simplify the calculations in MOND. 



12] and recent [l9( physics literature formulated Milgrom's 



(22) 



However, this formulation is incorrect except for a sharp-break function pL = min(l, a/00), 
so a = max(ajv, A /aoOjv). At this particular situation, 



p, \a/a ) = a /a = y/a a N = fj, l/2 (a N /a ) 



(23) 



which in general does not hold for other choices of jl(x). 

In Fig. [U we compare | V$| of different forms. We also plot the results of the standard 
form and simple form of the formalism adopted in 13Ferreras et al. 19j. For the same form 
and mass, the MONDian gravity of ours is always stronger than that of 13Ferreras et al. 

ET 



C. Doubled Image systems 

In realistic strong lensing systems, light rays often penetrate the mass distribution of 
lens such that point mass model is not appropriate. Since MOND is supposed to have 
little dark matter, we adopt the Hernquist model of elliptical galaxies [26] for our lenses. 
The Newtonian gravity is given by |V$tv| — GM j '(g + rh) 2 ■ In the following calculations we 



assume r^ = 0.551r e , where r e is the effective radius from surface brightness observation 26 ]. 

As shown in last subsection, seeking solution a = V$ of Eq. depends on the choice 
of p. In general we can express |V$| = g(V$Ar), and the deflection angle produced by a 
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FIG. 1. The strength of the MONDian gravity in different forms of £l(x). 
spherical lens (Eq. f JT2]) ) can be written as 

4GM, 

A<p = /, (24) 

where 

f = L GM^-4Y" dQ - (25) 

Moreover, since what we can measure in sky are not deflection angles but positions of the 



projected images, it is useful to define 6 = Qq/Dl and Be = V 'AGMDlsi 1 DlD$ , where 
Dl,D$ and DlS are distances from the observer to lens, observer to source and lens to 
source, respectively. Here $e is called the Einstein radius. For a spherical strong lens, two 
images are located on both sides of the lens (9 + and 6L). The corresponding lens equations 

8 



are [cf., |l3] 



D 



8% f 



s 



s v- 
where is the source position. These two lens equations can be combined into 

0+040+ + 0_) 

E (04_ + 044- 

Therefore, with the observed positions of the two images, + and 0_, we are able to infer 
the total mass of the lens by computing 9e ■ 



(26) 
(27) 

(28) 



III. RESULTS 

We apply our formalism to 10 double-image lenses from the CASTLES Catalogue [27J, 

of which the stellar masses have been estimated by stellar population synthesis model with 

different initial mass functions (IMFs) [2l|. Table 1 lists our result in conventional ACDM 

cosmology (fl m = 0.3, f^A = 0.7, h = 0.7). The lensing mass Ml estimated from the 

Bekenstein's form is the smallest, then the simple form and then the standard form, on top 

□ 

of them is the standard form by 13Ferreras et al. [19J. Furthermore, comparison of Ml 
from GR and TeVeS supports the idea that the mass discrepancy between the Newtonian 
and the MONDian paradigm can be quite significant. 

For a more consistent analysis, lensing in TeVeS should be studied in the MONDian 
cosmology - z/HDM (Sl b = 0.05, Vt u = 0.17 fi A = 0.78, h = 0.7) (tJ. We list the total mass of 
the 10 systems for the three forms of ft(x) in the bracket of columns 3 to 5 (in bracket) of 
Table 2. It appears that in TeVeS, lensing mass in the z/HDM cosmology is about 5% less 
massive than in ACDM. 



12Ferreras et al. [21[ estimated the aperture stellar mass by two IMFs: 03Chabrier [28] 
and 23Salpeter [29], and the results are listed in columns 6 & 7 of Table 2, respectively. For 
comparison, we compute the mass enclosed inside the truncated radius given in 12Ferreras 



et al. [2l|] from the lens total mass. The result is listed in columns 3 to 5. 

When comparing with the stellar mass (Ms a i p ) from Salpeter's IMF, we find that, except 
for BRI0952 — 0115 (where lensing gives smaller mass), all masses derived from the simple 
form are within the uncertainty of Mg a i p . However, when comparing with stellar mass 



TABLE I. Total mass of lenses (1O 1O M ) in ACDM 

TeVeS GR 



Lens 


Bekenstein 


Simple 


Standard 


FSY08.f. 




FSY08 


Q0142 - 100 


19.39 


24.32 


28.34 


29.28 


32.29 


32.37 


HS0818 + 1227 


29.59 


37.64 


44.68 


46.31 


50.80 


50.99 


FBQ0951 + 2635 


2.20 


2.72 


3.04 


3.82 


3.30 


4.07 


BRI0952 -0115 


2.59 


3.35 


4.05 


6.62 


4.60 


7.33 


Q1017 - 207 


6.41 


8.17 


9.69 


9.04 


10.95 


9.93 


HE1104 - 1805 


64.53 


82.99 


99.76 


103.17 


112.67 


112.93 


LBQ1009 - 025 


11.65 


14.65 


17.02 




19.28 




B1030 + 071 


17.55 


21.54 


24.40 




27.08 




SBS1520 + 530 


17.92 


22.68 


26.51 




29.57 




HE2149 - 274 


14.33 


18.14 


21.36 




24.14 





a Ferreras et al.(2008) 

(Mchab) from Chabrier IMF, Bekenstein's form is a better choice: 8 objects agree with 

Figs. |2]&|3]show the mass ratio of TeVeS to GR for the 10 lensing systems. For comparison 
of mass differences between the three forms, all ratios are plotted against |V$|/oo of the 
simple form. We found the ratio increases slightly with |V$|/oo. This does make sense 
because a smaller |V$|/oo means the system is closer to the MOND regime, and a larger 
mass discrepancy is expected. 

IV. DISCUSSION 

In our analysis, the simple form of £l(x) (i.e., a = 1, rj = 1 in the canonical form (1191)) 
yields the most reasonable lensing mass in TeVeS: in 9 of 10 systems the lens mass Ml is 
within the uncertainty of Ms a i p , the mass estimated from population syntthesis model with 
Salpeter's IMF. When compare with result estimated from Chabrier's IMF, 4 systems are 

10 



TABLE II. Aperture mass (and total mass) of lenses (1O 1O M0) in z/HDM 

TeVeS IMF (FSW05.D 



Lens 


|V$ s |/a 


Bekenstein 


Simple 


Standard 


Chabrier 


Salpeter 


Q0142 - 100 


5.85 


10.79 (18.36) 


13.66 (23.20) 


16.05(27.31) 


20.9?°;* 


18.^X3.2 


HS0818 + 1227 


5.67 


18.14 (28.30) 


23.39 (36.17) 


27.79 (43.35) 


16.2&g 


20.8ii 


FBQ0951 + 2635 


10.23 


1.54 (2.16) 


1.91 (2.67) 


2.15 (3.01) 


L - L 0.5 




BRI0952 -0115 


5.28 


2.01 (2.48) 


2.59 (3.21) 


3.19 (3.93) 


3.5^7 




Q1017 - 207 


5.16 


2.45 (5.89) 


3.22 (7.55) 


3.81 (9.15) 


4.3} 3 4 ° 


6.4^o 


HE1104 - 1805 


4.98 


45.17 (59.58) 


58.44 (77.08) 


71.78 (94.68) 


on o51.2 
zz -°12.7 


36.6l;I 


LBQ1009 - 025 


5.40 


7.71 (10.79) 


9.76 (13.67) 


11.53 (16.15) 


c r.7.9 
d.0 4 2 




B1030 + 071 


8.06 


9.76 (16.61) 


12.06 (20.51) 


13.80 (23.47) 


10.6& 3 


14.5^ 


SBS1520 + 530 


6.03 


11.91 (16.67) 


15.20 (21.28) 


18.08 (25.31) 






HE2149 - 274 


5.86 


7.04 (13.58) 


8.96 (17.28) 


10.67 (20.58) 


4-6| 6 7 


6-91:8 



Ferreras et al.(2005) 



outside the uncertainty of Mchab, but the mass difference is within 17~^q^%. Moreover, the 
simple form yields an average of 28% mass discrepancy between the inferred lensing mass of 



TeVeS and GR. This is close to the value 30% estimated by 24Sanders [30j. For comparison, 
Bekenstein's form gives an average of 44%, and provides a better fit to lower-mass IMFs. 
This consists with the conclusion of 20Napolitano et al. {31 1. 

We have to keep in mind that the uncertainty in estimating the mass of elliptical galaxies 
from IMF is still quite large. It is far from mature to claim that the MONDian paradigm 
favors which choice of p-(x) in elliptical galaxies. However, our study shows that like spiral 
galaxies, elliptical galaxies can be used to distinguish different forms of ft(x) as well. We 
believe that along side with an independently precise measurement of mass of elliptical 
galaxies is available, strong gravitational lensing can offer us a window to study the form of 
fi(x). We explicitly write down a formalism for this application in the future. 

We conclude that these systems do not show any apparent need of ad hoc dark matter in 
elliptical galaxies. On the other hand, with a simple spherically symmetric lens model, TeVeS 
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FIG. 2. Mass differences between Mqr and MpeVeS in three different forms of fi(x). Estimated 



stellar mass using Salpeter's IMFs 



21 



is shown for comparison. 



is able to fit these system reasonably well. This differs from the conclusion of 13Ferreras et al. 

n n n 

[19] but sides with 34Zhao et al. [Uj and 26Shan et al. [16] , where a constant mass-to-light 
ratio is assumed. 

In 13Ferreras et al. jl9|. they concluded that MOND might have problems in explaining 
galactic lensing system. Ours analysis differs from their conclusion. The reason could be 
that they have applied the formalism of Eq. (122|) in their calculation. Indeed, this might also 
explain why they have found a considerable discrepancy between MOND and TeVeS 32 ]. 
In 16Mavromatos et al. 32] and 14Ferreras et al. 20], the full relativistic equations with 
vector were solved for the first time. They discussed the result of two forms of fi proposed 



by 01 Angus et al. 



23] . Their choices of \x are identical to the simple form and Bekenstein's 
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FIG. 3. Same as Fig. [21 except that the estimated stellar mass used Chabrier's IMFs 21]. 



form in this paper, but a gap still exists between our result and 16Mavromatos et al. 32 ]. 
The disparity might be due to the fact that we choose Hernquist model rather than NFW 
model for the lens. 

Our analysis also shows that the simple form seems to give the most reasonable angle 
of deflection in strong lensing, which coincides with the studies on the dynamics of spiral 
galaxies jsl, [22I 24]. It is quite reasonable because the deviation of TeVeS from GR in strong 
lensing is only due to the change of gravity (Possion equation). So we do expect a consistent 
conclusion between gravitational lensing and dynamical analysis. 

We should point out that our analysis is based heavily on the assumption that the lenses 
in these systems are spherical and quasi-static, so that we can use Eq.(J2J) to derive the 
lensing equation. However, the assumption of spherical distribution is obviously insufficient 
for most known lenses. Moreover, the choice of mass model for lens may also affect the 
inference of the total mass of the lens. Thus, it is worthwhile to study how the deviation 
from this assumption would affect our conclusion here. We note that while a spherical 



13 



Hernquist model is a reasonable assumption for some lenses, it is always a poor assumption 
for modeling rotation curves of spiral galaxies which have often exponential profiles in their 
disks, and the gravity is enhanced in the disk plane. In this regard the recent claim of 
inconsistency of lensing and galaxy rotation curve in MOND [20j has yet to be corrected for 
these systematic effects. Systems in cluster environment which are not yet in full d yna mical 



equilibrium it might be important to include non-trivial effects of the vector field 33|, |34|. 
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